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Abstract 

Scalar form factor of the nucleon is calculated in the framework of light cone 
QCD sum rules, using the most general form of the baryon current. Using the result 
on scalar form factor of the nucleon, the nucleon-scalar a and ao meson coupling 
constants are estimated. Our results on these couplings are in good agreement with 
the prediction of the external field QCD sum rules method. 
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1 Introduction 



In spite of the great success of QCD theory of strong interactions, the quark structure 
of scalar mesons is still an open question and widely discussed in literature (see pQ for a 
review) . There are many scenarios for the scalar meson structures like two-quark qq, four- 
quark qqqq states, meson-meson bound states and even a scalar glueball. It is very likely 
that the scalar meson structure is not made of such simple structures, but are superpositions 
of these contents. Therefore different scenarios can give quite different predictions on the 
production and decays of scalar mesons which can be tested in experiments. One efficient 
way to make progress about the nature scalar bosons is to study their role in the various 
two-baryon reactions (NN, AN, EN, AA, EE). 

Our aim in this paper is to calculate the scalar form factor of the proton in the framework 
of light cone sum rules (LCSR) [2l [3j using the most general form of the proton current 
and then relate it to the nucleon-scalar a or a® mesons coupling constants. This method 
is an alternative approach to the traditional sum rule ^j. LCSR is based on the operator 
product expansion (OPE) near the light cone, which is an expansion of the time ordered 
product over the twist rather than the dimension of the operators. Main contribution of 
this approach comes from the operators having lower twist. The main ingredients of the 
LCSR are the wave functions of hadrons which define the matrix elements of non-local 
operators between the vacuum and one hadron states. 

This method applied to study the nucleon electromagnetic form factors [5], and the 
weak A& — > p£9g decay [6] . The higher twist wave functions of baryons have been obtained 
recently in [TJ. 

The plan of this article is as follows: In section 2 we derive the sum rule for the proton 
form factor due to the scalar current. Section 3 is devoted to the numerical analysis and 
discussion. 

2 Sum rules for the proton form factor due to the 
scalar current 

In the present section we derive sum rules for the scalar form factor of the proton using the 
most general form of nucleon interpolating current. For this aim we start with the following 
correlation function: 



where J(x) = uu + (— 1) dd is the interpolating current for the scalar a (J = 0) and ao 
(1 = 1) mesons, respectively, and r] is a suitable interpolating current with nucleon quantum 
number. Few words about the interpolating current of the nucleon are in order. There exist 
the following currents with the isotopic spin 1/2, without derivative terms (see [8]): 




(1) 



2 




(2) 



abc 



[ (u Ta (x)C^u\x)) i4d c (x) - (u Ta (x)C^d\x)) i4u c (x)] , (3) 



1 



where A\ = I, A\ = 75, A\ = 75, A\ = j3, and in C is the charge conjugation operator, 
a, b, c are the color indices, and z in Eq. ([3]) is a light-like vector with z 2 = 0. In Eq. (j2J), 
(3 = —1 corresponds to the Ioffe current [9]. Note that in [5] 772(0;) is modified as 

V3 (x) = e abc [(u Ta {x)C^u\x)) ^d c (x) . 

The nucleon scalar form factor which is obtained using the current 773(2:) in the light cone 
QCD sum rule approach is studied in [10]. However, the current 773(2;) couples also to 
the states with isospin 3/2, and these contributions should be eliminated from the final 
result, which is not done in [ID] . Moreover, terms that are proportional to and tfjt in 
the coordinate space which give contribution to the considered structure in [H] . are also 
neglected. In the present work 771 is chosen as the general form of the interpolating current. 

Let us firstly calculate the hadronic part of the correlator function. By inserting the 
complete set of states between the currents in Eq. ([I]) with quantum numbers of the 
corresponding nucleon, we get the following hadronic representation of the correlator 

Wp,q)- <%W)> W)!- 7 !^)) 1 y (0\vW))W)\J\N{p)) ^ 
m 2 N — p' 2 ^ ml — p' 2 

where p' = p — q, and q is the momentum carried by the scalar current. The second term in 
Eq. (j3J) takes into account higher states and continuum contributions and hadrons h form 
a complete set of baryons having the same quantum number as the ground state baryon B 
(here B represents p and n). The matrix elements in Eq. (j3J) are defined as 

(0\V\N( P ')) = \ N u N (p>) , (5) 
(N(p')\J\N(p)) =g(Q 2 )u(p')u(p) . (6) 

Using Eqs. (SJ) — (JSJ) and performing summation over nucleon spin, for the phenomenological 
part of the correlator we obtain 

n(p,g)= A ^ (g ; )( ^; +m " Mp) + -- - , (7) 

m N - VP - Q) 

where • • • represents contributions from the higher states and the continuum. It follows from 
Eq. (J7|) that the correlator function contains three different structures, namely structures 
~A ~/4 an d structure proportional to the unit operator /. In further analysis we will 
choose the structure 

Here we would like to make the following remark: The nucleon-scalar meson S coupling 
constant gNNS can be related to the scalar form factor g(Q 2 ), and for this aim we define 
the following coupling constant 

(N(p')S\N{p)) = g NNS u(p')u(p) . (8) 

It follows then from Eq. ([I]) that 

u( , _ (0\J\S) (N(p')S\N(p))(0\ V \N(p')) 
[P,q> (m 2 N - p' 2 )(m 2 - q 2 ) ' 1 ' 
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The matrix element of the scalar current J between the vacuum and the scalar meson state 
is defined as 



(0|J| S) = \ s 



(10) 



Using Eqs. (!H!)- (|T0l) and comparing with Eq. (J7J), we get the following relation between 
nucleon scalar form factor and Qnns 



g{Q) 



2 . 9nns^s 



m 2 s + Q 2 



At large Euclidean momenta p' 2 = (p — q) 2 and q 2 = —Q 2 , the QCD part of the correlator 
can be calculated from the explicit form of interpolating currents J{x) and 77(0) , from which 
we obtain 



n = ~ J d^xe iqx J2 1 {ca[ 



cry 



AiS(-x) 



Pf3 



Ae abc (Q 



T{<(0)^(x)^(0)} 



N 



+ (Ai) [cA{S(-xj 



—X) 



7/3 



Ae abc ( 



Ae abc ( 



CV7 



T{<(0)4(a;)^(0)} 
t{<(0)^(0)}^( x ) 



N 




(12) 



where S(x) is the light quark propagator. It follows from Eq. (TT2|) that for calculation of 
the correlator from QCD side we need to know the explicit form light quark propagator 
and the nucleon distribution amplitudes (DAs). The full propagator of the massless quark 



is 



S{x) 



i yt (qq) 



2n 2 x 4 



12 



1 + 



2 2 ' 

VTLqX 

16 



dv 



\Qtx 2 x a 



4ir 2 x 2 



■ (13) 



The contribution ~ G^ u can give rise to four and five particle nucleon DAs. These DAs are 
usually small [T2"l [T5] and can be neglected. 

The second term in Eq. ( TIB"]) does not give any contribution to the sum rule after the 
Borel transformation. So, only the first term is relevant for our calculations. 

As has already been noted, in order to calculate the non-perturbative contributions to 
the theoretical part of the correlator function one needs to know the matrix element of 
e abc (0\Ua(aix)u b 3 (a2x)d'^(a 3 x)\N(p)) between the proton and vacuum states. This matrix 
element of the non-local operator is parametrized in terms of the nucleon DAs [H Q2J [EU 
[T5l [T6| [T7] . Their explicit forms are given in the Appendix-A. 

Taking into account the three valence quark light-cone DAs up to twist-6 and performing 
integration over x in the coordinate space, we finally obtain the following expression for the 
correlator function: 



n 



1 

-m N 



dtn 



l ~ t2 dtv 4- fa 



4 " Jo ~~" Jo (q - pt 2 ) 2 



{2(1 + (3)V 2 (U) + 2[3 T 1 + (3 ± l)P]Vs(ti) 



±2(l + /?)P 1 (t,)-2[2±l 
+ 2[3±1 + (3t I) 13} A 3 (U] 



\-[2 T 1 - (2 ± l)(3]T 2 (U) + 12(1 



/3)T 7 (U 
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+ [6Tl-(6±l)/TO*i)} 



m 



A' 



dr / <iA 



l-ta 



1— PT 



(q - pr) 4 



{±(1-/3)^) 



T (1 - + 3[2 =f 1 — (2 ± l)/3]T 8 (t t )} 



dA / dU 



i-t 2 



dtr 



:{[2 T l + (2±l)P]V 2 (t l 



(q-p\) 2 

+ [2 ± 1 + (2 t l)0\MU) + [2 T 1 - (2 ± 

±2(1 + /3)T 4 (^)} (14) 
where any function F(ti) means F(tj) = F(ti,t 2 , 1— ti, t 2 ), and functions with ~ are defined 



as 



v 2 {u 

A 2 (t i 
A 6 {U 

%(ti 



V 1 (t i )-V 2 {t i )-V 3 {t i ) , 
-Aifc) + A 2 (U) - MU) , 

A^U) - A 2 (U) + A 3 (U) + A A {U) - A 5 (U) + A 6 (U) , 
T 1 (t i )+T 2 (t i )-2T 3 (t i ) , 
T 1 (t i )-T 2 (t i )-2T 7 (t i ) , 

2T 2 {U) - 2T 3 (U) - 2T A {U) + 2T 5 (^) + 2T 7 (U) + 2T 8 (t;) , 
-TxiU) + T 2 (U) + T 5 (ti) - T 6 (ti) + 2T 7 (U) + 2T 8 (U) . 

Explicit expressions of A a (ti), T a (ti) and V a (ti) can be found in [17]. Note also that in 
deriving Eq. ( fT4l) we omit the terms proportional with the unit operators which do not give 
contribution to the structure Upper (lower) sign in Eq. ( fill) corresponds to the scalar 
form factor due to the scalar current with the quantum numbers of a (ao). 

In order to suppress the continuum and higher state contributions we need to apply Borel 
transformation to the phenomenological and theoretical parts of the correlation function 
with respect to the variable (q — p) 2 . 

For the theoretical part the Borel transformation and the continuum subtraction can 
be done by using the following substitution rules 



dx 
dx 



p(x) 



(q — xp) 2 
(q — xp) 4 



X 



1 

M 2 



dx PM e s/M> + 



Q 2 + x\M 2 



,-so/M 2 



(15) 



where 



1 -x)M 2 + 



x 



X 



Q 2 



(Q 2 + s - m 2 N f + 4m^Q 2 - (Q 2 + s - m 



N) 



and Q 2 = 

g(Q 2 ) 



2mjv62 

-q 2 . Finally, we get the following sum rule for the scalar form factor. 



e m%/M* r 



A 



N 



-m N 



dU 



-s/M' 



{2(1 + /3)V 2 (U) + 2[3 T 1 + (3 ± l)/3]V 3 (ti 
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± 2(1 + /3)Px(U) - 2 [2 ± 1 - (2 t iWxiU) + 2(1 + (3)A 2 (U) 

+ 2 [3 ± 1 + (3 T l)/3]^ 3 (^) + [2 T 1 - (2 ± l)/3]T 2 fo) + 12(1 - /3)T 7 fo 

+ [6Tl-(6±l)/3]^(ti)} 



+ 



A' 



£dr£d\£dt2 / 1 "* 2 ^i^ e " /M2 { ± (l-/3)V 6 (*i) 



T (1 - /3)A(*i) + 3 [2 T 1 - (2 ± l)/3]T 8 (t,)} 



+ 



Q 2 + xlm 2 N Ji 



■t'o 



<iA / <it 2 



i-t 2 



tM 2 



s / M2 {±(l-/3)V 6 (t.) 



=F (1 — + 3 [2 T 1 - (2 ± l)/3]T 8 (t t )} 



i dU r 1 -** 



^ / dt 1 e^ M2 {[2 T l + (2±l)P}V 2 (t l 

x t 2 J L 

+ [2 ± 1 + (2 t l)0\MU) + [2 T 1 - (2 ± l)/3]T 2 (t 4 ) 
±2(l + /?)T 4 (t i )} 



(16) 



It follows from Eq. ( fl6l) that for determining the scalar form factor, the residue A at of the 
nucleon needs to be known. The residue A at is obtained from the mass sum rules for the 
baryons (see for example [IB] ) 



X 2 Ne -m%/MZ 



-£ 2 (a;)(5 + 2/3 + /3 2 ) + 5 ' 



256tt 4 " 
2AM 2 



6 



6(l-/? 2 )(<M) + (-l+ / 3 2 )(?m) 



(tZtt) -12(1- (3 2 ){dd) + {-l + (3 2 )(uu) 



(17) 



where x = sq/M 2 , and 



£ n (x) = l-e^ 



n k 



k=0 



which describes the subtraction of higher states and continuum contributions. 



3 Numerical Analysis 

In this section we present the numerical analysis for the scalar form factor g(Q 2 ). It 
follows from Eq. ffl6|) that the main ingredients of LCSR are the DAs, which involve 
non-perturbative parameters fjy, Aj, A 2 , ff, / 2 , Af, Vf (see Appendix-B). 

In the numerical calculations we use three different sets of the DAs (for more details, 
see [5] and [17]). 

1) QCD sum rules based DAs, 

2) Asymptotic DAs, 

3) A model for the nucleon DAs. The parameters in this case are chosen in such a way 
that, the form of DAs describes very well the existing experimental data on nucleon form 
factors. 
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The values of the above-mentioned non-perturbative parameters for these cases at fi = 
1 GeV scale are given in [5] and [17]. Our numerical calculations show that the values of 
the scalar form factors of nucleon due to the scalar current in QCD based DAs case are 
larger compared to that of the other two scenarios, and, for customary reasons we present 
our numerical results only for this case. 

In the problem at hand we have three auxiliary parameters, namely, Borel mass square 
M 2 , continuum threshold Sq and /3 in the interpolating current of nucleon. Obviously, any 
physically measurable must be independent of the auxiliary parameters. So, in the first 
hand we have to find the appropriate regions of these parameters, in where the form factor 
is independent of them. The mass sum rules analysis shows that optimal value for the 
continuum thresholds is in the region sq(2 4- 2.5) GeV 2 . 

For this aim we consider the following three-step procedure. In the first step we try 
to find the working region of M 2 where scalar form factor is independent of M 2 at fixed 
values of s and (3. In Figs. (1), (2), (3) and (4) (Figs. (5), (6), (7) and (8)) we present 
the dependence of the scalar form factor induced by I = (/ = 1) scalar current on M 2 at 
different fixed values of Q 2 and [3 at s = 2.0 GeV 2 and s = 2.5 GeV 2 , respectively. 

From these figures we see that, in both cases, the scalar form factor is practically inde- 
pendent of M 2 at different values of Q 2 , (3 and Sq. To decide on an upper bound for M 2 we 
require the contribution of continuum be less than the contribution of continuum-subtracted 
sum rules and the lower bound can be determined by requiring that the contribution of the 
highest power of 1/M 2 to be less than, say, 30% of the higher powers of M 2 . As a result 
of our analysis we find that the region 1 GeV 2 < M 2 < 2.5 GeV 2 satisfies both conditions 
simultaneously. 

The parameter (3 is an auxiliary parameter, and hence, the physical quantities should be 
independent of it. For this reason we need to find the working region of (3 where the scalar 
form factor is independent of it. In determining this region of f3 two conditions should be 
satisfied: Eq. (ITT)) and its logarithmic derivative with respect to 1/M 2 should be positive 
and independent of (3. From an analysis of Eq. (ITT)) it is found that —0.5 < (3 < 0.5 is 
unphysical [18J. It is further shown in [18] that the working region of f3 where the second 
condition is satisfied is obtained for —0.6 < cos# < 0.3, where tan# = (3, which yields 
(3 < -1.3 and (3 > 3.3. 

We perform calculations depicting the dependence of the scalar form factors on cos 9 at 
fixed values of s (s„ = 2.0 GeV 2 and s = 2.5 GeV 2 ), M 2 (M 2 = 1.0 GeV 2 and M 2 = 
1.5 GeV 2 ) at three fixed values of Q 2 {Q 2 = 3.0 GeV 2 , Q 2 = 4.0 GeV 2 and Q 2 = 5 GeV 2 ); 
and confirmed that in the regions when (3 < —1.3 and (3 > 3.3 ,the form factor is practically 
independent of (3. So, the common region where mass sum and form factor are independent 
of (3 is the region —0.5 < cos 9 < 0.3, which corresponds to (3 < —1.3 and (3 > 3.3 and which 
we will use in further discussions. Note that, the Monte Carlo analysis that is performed in 
[T9] predicts the optimal value of (3 to be (3 = —1.2. Our lower bound is slightly different 
compared to this value. 

We perform the OPE in the light cone with large Euclidean Q 2 and (p — q) 2 momenta, 
where the scalar from factor g(Q 2 ) can reliably be determined at the range Q 2 > 2 GeV 2 . 
In the region Q 2 < 2 GeV 2 OPE becomes questionable and for this reason in determining 
g(Q 2 ) we restrict ourselves to the region Q 2 > 2 GeV 2 . 

In Figs. (9), (10), (11) and (12) we present the dependence of g{Q 2 ) induced by / = 
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0(7 = 1) current on Q 2 at fixed values of (3 which lies in the above-mentioned region of j3; 
at M 2 = 1 GeV 2 and s = 2 GeV 2 , s = 2.5 GeV 2 , respectively, for the central values all 
parameters. We observe that the value of g(Q 2 ) seems to be insensitive to the choice of so- 
We further see from these figures that for positive (negative) values j3, the sign of the form 
factor is positive (negative). Therefore determination of the sign of g{Q 2 ) can give valuable 
information about parameter f3. 

Having the obtained the Q 2 dependence of the form factor, we can determine Qnns- It 
follows from Eq. (flOl) that 

(m 2 s + Q 2 )g(Q 2 ) 

9nns — 7 • 

*s 

The residues \ a and A ao of the scalar currents are determined from two-point QCD 
sum rules analysis in [201 [2T], which is predicted to have the value A^ = 0.2 and A ao = 0.3. 
Using the results of Figs. (3) and (4), for the scalar form factor gNNa we get 

6 ± 1 for /3 = 5 , 



\dNNcr\ 



\9NNa \ 



3 ± 2 for p = -5 ; P = -3 
5.5 ± 1.5 for p = 5 , 



7 ± 1 for p = -5 ; p = -3 . 

For completeness we present the results for gNNS coupling for the Ioffe current case, 
despite that P — — 1 lies outside the working region of p. Our results are 

\9NNa\ = 12 ± 2 , 
\9NNa \ = 11 ± 2 . 

Note that his coupling constant in framework of the external-field QCD sum rules 
method using the Ioffe current is evaluated in [22], and it is obtained that gNNa = 14. 4 ±3. 7. 
This same coupling constant is also shown to have the value gNNa — 16.9 pQ from a fit to 
the NN scattering data. 

A comparison of our result for p = — 1 shows that, within the error limits, there is a 
good agreement between our result and the above-mentioned one. 

In conclusion, we have calculated the nucleon scalar form factor in the framework of 
light cone QCD sum rules, using the most general form of the baryon current. Using 
the obtained result on scalar form factor, we estimate gNNS- Our result in regard to the 
coupling constant gNNa is in a good agreement with the prediction of the external field 
QCD sum rules method for the Ioffe current. Measurement of the nucleon scalar meson 
coupling constant can give valuable information about the quark structure of scalar meson 
as well as about the value of parameter p. 
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Appendices 



A Nucleon distribution amplitudes 

In this appendix we present the three-quark distribution DAs from twist-3 to twist— 6 (for 
more details, see [Tj and [T7]). 

The matrix element of the non-local operator e abc (0\u^ l (aix)u b /3 (a2x)d'^(a3x)\N(p)) be- 
tween the proton and vacuum states is parametrized in terms of the nucleon DAs as follows: 



4e afec (0 



u a a (aix)u b p(a2x)d c 1 (a^x) 



P 



S 1 m N C a ^('y 5 N) 1 + S 2 m 2 N C af 3(fy 5 N), 



+ V 1 m N ( l5 C) a ^N 1 + V 2 m 2 N { lb C) aP {^N) 1 + Vi + 



2 2 
X m N^iM 



+ V 2 m N {PC) af3 {fr 5 N) 7 + V 3 m 7 v(7 At CV(7^7 5 A0 7 + V 4 m 2 N {^C) a p{ lb N)^ 
+ V B m%('y lt C) a f,(i<T fa 'x v 'y B N)i + V 6 m^(^) a/3 (^7 5 Ar) 7 

+ (At + J (P 75 C) a/3 iV 7 + A 2 m N (P l5 C) af3 (^N\ + A 3 m N ( lfll5 C) a ^N), 

+ A 4 m 2 N (^ 5 C) a pN^ + A b m 2 N {y^r > C) aP {ia lxv x v N) 1 + A 6 m 3 N (fy 5 C) al3 {^N) y 
+ ( Tl + ^3irA (P^CUiYlsN), + T.m^P^a^CUi^N), 



4 



+ Tgm^cvCVK^s A0 7 + T A m N (P' / a^C) a(3 (cr^x p -f 5 N)^ 
+ T 5 m 2 N {x v ia flu C) aP {^ 5 N) ry + T 6 m 2 N (x»P u ia^C) a(3 (fy 5 Ny 
+ r 7 m 2 N (a^C) af3 (cr^ ^ 5 N\ + %m 3 N (x v a^C) af3 ((j'" , x pl5 N) 



7 ■ 



(A.1) 



The calligraphic DAs do not have definite twist and can be related to the one with definite 



twist as: 














<Si — Si , 






(2P 


x] 


s 2 = 


Si — s 2 , 


Vi = Pi, 






(2P 


x] 


v 2 = 


Pi - Pi , 


Vi = Vi, 






(2P 


x] 


V 2 = 


Vi-V 2 -V 3 , 


2V 3 = V 3 , 






(4P 


x] 


V 4 = 


-2Vi + V 3 + V 4 + 2V 5 , 


(4P-x)V 5 = 




-v 3 , 


(2P 


x) 


2 v 6 = 


= -v x + v 2 + v 3 + V A + v 5 - v 6 , 


A 1 = Ai, 






(2P 


x) 


A 2 = 


-A 1 + A 2 - A 3 , 


2A 3 = A 3 , 






(4P 


x) 


Ai = 


-2A 1 - A 3 - A 4 + 2A 5 , 


(AP-x) A 5 = 


A 3 - 




(2P 


x) 


2 A 6 - 


= Ai - A 2 + A 3 + A 4 - A 5 + A 6 , 


t 1 = t 1 , 






(2P 


x) 


T 2 = 


Ti + T 2 - 2T 3 , 


2T 3 = T 7 , 






(2P 


x) 


% = 


Ti-T 2 - 2T 7 , 


(2P-x)T 5 = 


-Ti 


+ T 5 + 2T 8 , 


(2P 


x) 


2 T - 


-- 2T 2 - 2T 3 - 2T 4 + 2T 5 + 2T 7 + 2T t 


(AP-x) T 7 = 


T 7 - 


T 8 , 


(2P 


x) 


2 T - 
1 8 - 


-- -Ti + T 2 + T 5 - T 6 + 2T 7 + 2T 8 . 



These DAs can be written as 

F(cLi p ■ x) = J dxi dx 2 dx 3 5{x x + x 2 + x 3 - 1) e - ip - xS ^ a * F( Xi ) . (A.2) 

These light cone DAs are scale dependent and can be expanded next-to-leading spin 
accuracy with the conformal operators whose explicit expressions are, 

Vi(xi,/i) = I20xix 2 x 3 [(f)l(fi) + - 3x 3 )} , 
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V 2 (xi,fi) = 24x1x2 [02 (A*) + 04 (A*)(l -5ar 3 )] , 
V^x,,^) = 12z 3 {^(aO(1 -x 3 ) + {n)[x\ + x 2 -x 3 (l -x 3 )] 

+ ipi(/J>)(l - x 3 - 10xix 2 )} , 
Vi(Xi,fl) = 3{^g(^)(l -x 3 ) + ip^(fi)[2x 1 x 2 -0:3(1 -x 3 )] 

+ ^t^)[l-xz-2(x\ + xl)}} , 

V 5 (x t ,fi) = 6x 3 [0°(//) + 0+(ai)(1 - 2x 3 )] , 
Vs{ Xi ,n) = 2[0°(/x) + 0+(/i)(l - 3x 3 )] , 
Ai( Xi ,n) = 12Oxix 2 x 3 03 (/i)(x 2 — x±) , 
A 2 (xi,fj) = 24x x x 2 04 (aO(x 2 > 

A 3 (x,,/i) = 12x 3 (x 2 -x 1 ){[^°(/i)+^ 4 + (/i)]+ V 4 "(A*)(l-2x 3 )} , 

A 4 { Xi ,fi) = 3(x 2 -xi){ - ^°(a*) +"05"(a*)^3 + -0^(a*)(1 - 2x 3 )} , 
A 5 (^i,A*) = 6x 3 (x 2 -xi)05 (//) , 

A 6 (Xi,/i) = 2(x 2 - Xi)0g (//) , 

7\ (X,,//) = ^Ox^X^/l) + 1(03 - 3 + )(/i)(l - 3x 3 )] , 
T 2 (xi,n) = 24xix 2 [^(/i) + - 5x 3 )] , 

T 3 ( Xi ,fi) = Qx 3 {(£ + 0° + ^ 4 °)(Ai)(l - x 3 ) + (£ 4 - + 4 - i>l )(fi)[xl + x\- x 3 (l - x 3 )] 

+ (£ 4 + + <Pt + i>t) (A*)(l - ^3 - 10xix 2 )} , 
T 4 (x l? /i) = ^{(£ 5 ° + 0^ + V 5 °)(A*)(1 - *s) + fe" + 05 - ^)(/^)[2^iX 2 - x 3 (l - x 3 )] 

+ {it + 0+ + - x 3 - 2(x? + **)]} , 

Ts^,//) =6x 3 [e 5 °(/i)+e5 + (^)(l-2a: 3 )] , 
T 6 ( Xi ,ri = 2[0°(/x) + 1(0- - 6 + )( Ai )(l - 3x 3 )] , 

T 7 (xi,fi) = 6x 3 {(-£ 4 ° + <p° A + ^)(A*)(1 - a*) + (-C + 01 - ^4~)(A*)[^ + 4 ~ a*(l - a*)] 

+ ("^ + 01 + ^4 + )(/")(l - ^3 - lOXxX,)} , 

T 8 (x 4 , /x) = ^{(-£ 5 ° + 0° + ^ 5 °)(/i)(l - x 3 ) + (-£ 5 ~ + 5 - - V 5 -)(A*)[2x!X 2 - x 3 (l - x 3 )] 

+ + ^ + ^t) (A*)[l - - 2(x? + xl)]} , 
^(x,, //) = 6x 3 (x 2 - xx) [(f° + 0^ + ^ 4 ° + e 4 + + 0! + ^ 4 + )(/i) + (C + 01 - ^4 _ )(A*)(1 - 2x 3 ) 

S 2 (x,, /i) = ^(X 2 - xOt-^ + 4 + + + 05" - C)(/^3 

+ (e5 + + 05'+^5 + )(^)(l-2x 3 )], 

P 1 (x i , n) = 6x 3 (x 2 - Xl )[(e A - 02 - ^ 4 ° + C - 4 + - ^4 + )(a0 + (C - 04" + ^ 4 ")(a*)(1 " 2x 3 )] , 

P 2 (x l)Ai ) = ^(a; 2 -xO[^5 O + ^5 O -e5 O )(^) + (C-05 +^ 5 ")(^3 
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+ (£ 5 + -<^-^ + )0")(l-2z 3 )] , (A.3) 

where V\, A\ and T\ are leading twist-3; Si, Pi, V 2 , V 3 , A 2 , A 3 , T 2 , T 3 and T 7 are twist-4; 
S2, P2, V4, V5, A4, A5, T 4 , T 5 and T 8 are twist-5; and V 6 , A Q and T 6 are twist-6 DAs. 



B Non— perturbative parameters 

The coefficients <f%, (j)° 6 , <j>\, <f%, fg, ipl, ^5, #f, 03 , 04 , 04, V>4 » V»4 > £f , 64 > 
■0^, -05", ^5", 06 , 0g in the expansions of the DAs can be expressed in terms 
non-perturbative parameters /at, Ai, A2, ff, f£, Af, as follows (see [T7]): 
For the leading confirmal spin, 



<^5 > ^5 , 

of eight 



A? 



04 = 05 = 2 + Ai) 

$ = $ = 1*2, < = r 5 = \(f N ~xi) 



for the next-to-leading spin, for twist-3: 

3 " = ^f N Al, 

for twist-4: 



2^(1-3^ 



(B.l) 



(B.2) 



for twist-5: 



f N {3 -10^+^(3-10^ 



4 
5 r 



f N (l-2A^)-X l {l-2ff-Aff) 



f N {2 + f>A\-SV?) - Ai(2-5#-5tf) 



4 



16 



f N {2 -AI- 3Vi d ) - Ai(2 - 7/ x d + A" 

A 2 (4-15/ 2 d ), ^ = ^2(4-15^), 
lb 



5 
6 

5 r 
3 

5 r 

6 

5 

3 



/ w (3 + 4^)-A 1 (l-4/ 1 d ) 
f N (l-2AV-\i(f?-f?) j 
f N (5 + 2A u 1 -2Vf) - Ai(l-2#-2tf) 
f^-Al-W^ + X^fi-fi) , 



-A 2 (2-9/ 2 d ), e 5 " 



(B-3) 



(B.4) 
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and for twist-6: 



f N (l-4V*)-\ 1 (l-2f* 



/ JV (l+4A?) + A 1 (l-4/ 1 d -2/?; 



(B.5) 



The values of the non-perturbative parameters /at, Ai, A2, ff, f$i ^t-, can De found 
in [5] and [TT] . 
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Figure captions 

Fig. (1) The dependence of the scalar form factor of nucleon due to the scalar current 
with quantum numbers of a meson, on Borel parameter M 2 at So = 2 GeV 2 for three differ- 
ent values of Q 2 and f3: Q 2 = 3 GeV 2 , Q 2 = 4 GeV 2 , Q 2 = 5 GeV 2 and (3 = 1, = 3, (3 = 5. 

Fig. (2) The same as in Fig. (1), but at s = 2.5 GeV 2 . 

Fig. (3) The same as in Fig. (1), but for values of f3: (3 = — 1, — 3, — 5. 

Fig. (4) The same as in Fig. (2), but for values of (3: (3 = — 1, (3 = -3, (3 = -5. 

Fig. (5) The same as in Fig. (1), but for the scalar form factor induced by the scalar 
current with quantum number of a meson. 

Fig. (6) The same as in Fig. (5), but at s = 2.5 GeV 2 . 

Fig. (7) The same as in Fig. (5), but for values of (3: (3 = —1,(3 = —3, (3 = —5. 
Fig. (8) The same as in Fig. (7), but at s = 2.5 GeV 2 . 

Fig. (9) The dependence of the scalar form factor induced by the scalar current with 
a meson quantum numbers, on Q 2 at M 2 = 1 GeV 2 and s = 2 GeV 2 , at fixed values of 
P = -1.4, (3 = -3, (3 = -5, (3 = 5. 

Fig. (10) The same as in Fig. (9), but at s = 2.5 GeV 2 . 

Fig. (11) The dependence of the scalar form factor induced by the scalar current with 
a meson quantum numbers, on Q 2 at M 2 = 1 GeV 2 and s = 2 GeV 2 , at fixed values of 
P = -1.4, P = -3, P = -h, (3 = 5. 

Fig. (12) The same as in Fig. (11), but at s = 2.5 GeV 2 . 
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